The relationship between 11 different polygenic longevity scores, parental lifespan, and disease diagnosis in the UK Biobank.

Large-scale genome-wide association studies (GWAS) strongly suggest that most traits and diseases have a polygenic component. This observation has motivated the development of disease-specific "polygenic scores (PGS)" that are weighted sums of the effects of disease-associated variants identified from GWAS that correlate with an individual's likelihood of expressing a specific phenotype. Although most GWAS have been pursued on disease traits, leading to the creation of refined "Polygenic Risk Scores" (PRS) that quantify risk to diseases, many GWAS have also been pursued on extreme human longevity, general fitness, health span, and other health-positive traits. These GWAS have discovered many genetic variants seemingly protective from disease and are often different from disease-associated variants (i.e., they are not just alternative alleles at disease-associated loci) and suggest that many health-positive traits also have a polygenic basis. This observation has led to an interest in "polygenic longevity scores (PLS)" that quantify the "risk" or genetic predisposition of an individual towards health. We derived 11 different PLS from 4 different available GWAS on lifespan and then investigated the properties of these PLS using data from the UK Biobank (UKB). Tests of association between the PLS and population structure, parental lifespan, and several cancerous and non-cancerous diseases, including death from COVID-19, were performed. Based on the results of our analyses, we argue that PLS are made up of variants not only robustly associated with parental lifespan, but that also contribute to the genetic architecture of disease susceptibility, morbidity, and mortality.

The selection landscape and genetic legacy of ancient Eurasians.

The Holocene (beginning around 12,000 years ago) encompassed some of the most significant changes in human evolution, with far-reaching consequences for the dietary, physical and mental health of present-day populations. Using a dataset of more than 1,600 imputed ancient genomes1, we modelled the selection landscape during the transition from hunting and gathering, to farming and pastoralism across West Eurasia. We identify key selection signals related to metabolism, including that selection at the FADS cluster began earlier than previously reported and that selection near the LCT locus predates the emergence of the lactase persistence allele by thousands of years. We also find strong selection in the HLA region, possibly due to increased exposure to pathogens during the Bronze Age. Using ancient individuals to infer local ancestry tracts in over 400,000 samples from the UK Biobank, we identify widespread differences in the distribution of Mesolithic, Neolithic and Bronze Age ancestries across Eurasia. By calculating ancestry-specific polygenic risk scores, we show that height differences between Northern and Southern Europe are associated with differential Steppe ancestry, rather than selection, and that risk alleles for mood-related phenotypes are enriched for Neolithic farmer ancestry, whereas risk alleles for diabetes and Alzheimer's disease are enriched for Western hunter-gatherer ancestry. Our results indicate that ancient selection and migration were large contributors to the distribution of phenotypic diversity in present-day Europeans.

Polygenic profiles define aspects of clinical heterogeneity in attention deficit hyperactivity disorder.

Attention deficit hyperactivity disorder (ADHD) is a complex disorder that manifests variability in long-term outcomes and clinical presentations. The genetic contributions to such heterogeneity are not well understood. Here we show several genetic links to clinical heterogeneity in ADHD in a case-only study of 14,084 diagnosed individuals. First, we identify one genome-wide significant locus by comparing cases with ADHD and autism spectrum disorder (ASD) to cases with ADHD but not ASD. Second, we show that cases with ASD and ADHD, substance use disorder and ADHD, or first diagnosed with ADHD in adulthood have unique polygenic score (PGS) profiles that distinguish them from complementary case subgroups and controls. Finally, a PGS for an ASD diagnosis in ADHD cases predicted cognitive performance in an independent developmental cohort. Our approach uncovered evidence of genetic heterogeneity in ADHD, helping us to understand its etiology and providing a model for studies of other disorders.

Polygenic Analyses Show Important Differences Between Major Depressive Disorder Symptoms Measured Using Various Instruments.

BACKGROUND: Symptoms of major depressive disorder (MDD) are commonly assessed using self-rating instruments like the Patient Health Questionnaire-9 (PHQ-9) (current symptoms) and the Composite International Diagnostic Interview Short-Form (CIDI-SF) (worst-episode symptoms). We performed a systematic comparison between them for their genetic architecture and utility in investigating MDD heterogeneity. METHODS: Using data from the UK Biobank (n = 41,948-109,417), we assessed the single nucleotide polymorphism heritability and genetic correlation (rg) of both sets of MDD symptoms. We further compared their rg with non-MDD traits and used Mendelian randomization to assess whether either set of symptoms has more genetic sharing with non-MDD traits. We also assessed how specific each set of symptoms is to MDD using the metric polygenic risk score pleiotropy. Finally, we used genomic structural equation modeling to identify factors that explain the genetic covariance between each set of symptoms. RESULTS: Corresponding symptoms reported through the PHQ-9 and CIDI-SF have low to moderate genetic correlations (rg = 0.43-0.87), and this cannot be fully attributed to different severity thresholds or the use of a skip structure in the CIDI-SF. Both Mendelian randomization and polygenic risk score pleiotropy analyses showed that PHQ-9 symptoms are more associated with traits that reflect general dysphoria, whereas the skip structure in the CIDI-SF allows for the identification of heterogeneity among likely MDD cases. Finally, the 2 sets of symptoms showed different factor structures in genomic structural equation modeling, reflective of their genetic differences. CONCLUSIONS: MDD symptoms assessed using the PHQ-9 and CIDI-SF are not interchangeable; the former better indexes general dysphoria, while the latter is more informative about within-MDD heterogeneity.

Deep learning-based phenotype imputation on population-scale biobank data increases genetic discoveries.

Biobanks that collect deep phenotypic and genomic data across many individuals have emerged as a key resource in human genetics. However, phenotypes in biobanks are often missing across many individuals, limiting their utility. We propose AutoComplete, a deep learning-based imputation method to impute or 'fill-in' missing phenotypes in population-scale biobank datasets. When applied to collections of phenotypes measured across ~300,000 individuals from the UK Biobank, AutoComplete substantially improved imputation accuracy over existing methods. On three traits with notable amounts of missingness, we show that AutoComplete yields imputed phenotypes that are genetically similar to the originally observed phenotypes while increasing the effective sample size by about twofold on average. Further, genome-wide association analyses on the resulting imputed phenotypes led to a substantial increase in the number of associated loci. Our results demonstrate the utility of deep learning-based phenotype imputation to increase power for genetic discoveries in existing biobank datasets.

Phenotype integration improves power and preserves specificity in biobank-based genetic studies of major depressive disorder.

Biobanks often contain several phenotypes relevant to diseases such as major depressive disorder (MDD), with partly distinct genetic architectures. Researchers face complex tradeoffs between shallow (large sample size, low specificity/sensitivity) and deep (small sample size, high specificity/sensitivity) phenotypes, and the optimal choices are often unclear. Here we propose to integrate these phenotypes to combine the benefits of each. We use phenotype imputation to integrate information across hundreds of MDD-relevant phenotypes, which significantly increases genome-wide association study (GWAS) power and polygenic risk score (PRS) prediction accuracy of the deepest available MDD phenotype in UK Biobank, LifetimeMDD. We demonstrate that imputation preserves specificity in its genetic architecture using a novel PRS-based pleiotropy metric. We further find that integration via summary statistics also enhances GWAS power and PRS predictions, but can introduce nonspecific genetic effects depending on input. Our work provides a simple and scalable approach to improve genetic studies in large biobanks by integrating shallow and deep phenotypes.

Population-based Risk of Psychiatric Disorders Associated with Recurrent CNVs.

Recurrent copy number variants (rCNVs) are associated with increased risk of neuropsychiatric disorders but their pathogenic population-level impact is unknown. We provide population-based estimates of rCNV-associated risk of neuropsychiatric disorders for 34 rCNVs in the iPSYCH2015 case-cohort sample (n=120,247). Most observed significant increases in rCNV-associated risk for ADHD, autism or schizophrenia were moderate (HR:1.42-5.00), and risk estimates were highly correlated across these disorders, the most notable exception being high autism-associated risk with Prader-Willi/Angelman Syndrome duplications (HR=20.8). No rCNV was associated with significant increase in depression risk. Also, rCNV-associated risk was positively correlated with locus size and gene constraint. Comparison with published rCNV studies suggests that prevalence of some rCNVs is higher, and risk of psychiatric disorders lower, than previously estimated. In an era where genetics is increasingly being clinically applied, our results highlight the importance of population-based risk estimates for genetics-based predictions.

ADuLT: An efficient and robust time-to-event GWAS.

Proportional hazards models have been proposed to analyse time-to-event phenotypes in genome-wide association studies (GWAS). However, little is known about the ability of proportional hazards models to identify genetic associations under different generative models and when ascertainment is present. Here we propose the age-dependent liability threshold (ADuLT) model as an alternative to a Cox regression based GWAS, here represented by SPACox. We compare ADuLT, SPACox, and standard case-control GWAS in simulations under two generative models and with varying degrees of ascertainment as well as in the iPSYCH cohort. We find Cox regression GWAS to be underpowered when cases are strongly ascertained (cases are oversampled by a factor 5), regardless of the generative model used. ADuLT is robust to ascertainment in all simulated scenarios. Then, we analyse four psychiatric disorders in iPSYCH, ADHD, Autism, Depression, and Schizophrenia, with a strong case-ascertainment. Across these psychiatric disorders, ADuLT identifies 20 independent genome-wide significant associations, case-control GWAS finds 17, and SPACox finds 8, which is consistent with simulation results. As more genetic data are being linked to electronic health records, robust GWAS methods that can make use of age-of-onset information will help increase power in analyses for common health outcomes.

Multi-PGS enhances polygenic prediction by combining 937 polygenic scores.

The predictive performance of polygenic scores (PGS) is largely dependent on the number of samples available to train the PGS. Increasing the sample size for a specific phenotype is expensive and takes time, but this sample size can be effectively increased by using genetically correlated phenotypes. We propose a framework to generate multi-PGS from thousands of publicly available genome-wide association studies (GWAS) with no need to individually select the most relevant ones. In this study, the multi-PGS framework increases prediction accuracy over single PGS for all included psychiatric disorders and other available outcomes, with prediction R2 increases of up to 9-fold for attention-deficit/hyperactivity disorder compared to a single PGS. We also generate multi-PGS for phenotypes without an existing GWAS and for case-case predictions. We benchmark the multi-PGS framework against other methods and highlight its potential application to new emerging biobanks.

Circulating S100B levels at birth and risk of six major neuropsychiatric or neurological disorders: a two-sample Mendelian Randomization Study.

Circulating levels of the astrocytic marker S100B have been associated with risk of neuropsychiatric or neurological disorders. However, reported effects have been inconsistent, and no causal relations have yet been established. We applied two-sample Mendelian Randomization (MR) on the association statistics from genome-wide association studies (GWAS) for circulating S100B levels measured 5-7 days after birth (the iPSYCH sample) and in an older adult sample (mean age, 72.5 years; the Lothian sample), upon those derived from major depression disorder (MDD), schizophrenia (SCZ), bipolar disorder (BIP), autism spectral disorder (ASD), Alzheimer's disease (AD), and Parkinson's disease (PD). We studied the causal relations in the two S100B datasets for risk of these six neuropsychiatric disorders. MR suggested increased S100B levels 5-7 days after birth to causally increase the risk of MDD (OR = 1.014; 95%CI = 1.007-1.022; FDR-corrected p = 6.43×10-4). In older adults, MR suggested increased S100B levels to have a causal relation to the risk of BIP (OR = 1.075; 95%CI = 1.026-1.127; FDR-corrected p = 1.35×10-2). No significant causal relations were found for the other five disorders. We did not observe any evidence for reverse causality of these neuropsychiatric or neurological disorders on altered S100B levels. Sensitivity analyses using more stringent SNP-selection criteria and three alternative MR models suggested the results are robust. Altogether, our findings imply a small cause-effect relation for the previously reported associations of S100B and mood disorders. Such findings may provide a novel avenue for the diagnosis and management of disorders.

Polygenic analyses show important differences between MDD symptoms collected using PHQ9 and CIDI-SF.

Symptoms of Major Depressive Disorder (MDD) are commonly assessed using self-rating instruments like the Patient Health Questionnaire 9 (PHQ9, for current symptoms), and the Composite International Diagnostic Interview Short-Form (CIDI-SF, for lifetime worst-episode symptoms). Using data from the UKBiobank, we show that corresponding symptoms endorsed through PHQ9 and CIDI-SF have low to moderate genetic correlations (rG=0.43-0.87), and this cannot be fully attributed to different severity thresholds or the use of a skip-structure in CIDI-SF. Through a combination of Mendelian Randomization (MR) and polygenic prediction analyses, we find that PHQ9 symptoms are more associated with traits which reflect general dysphoria, while the skip-structure in CIDI-SF allows for the identification of heterogeneity among likely MDD cases. This has important implications on factor analyses performed on their respective genetic covariance matrices for the purpose of identification of genetic factors behind MDD symptom dimensions and heterogeneity.

Associations of psychiatric disorders with sex chromosome aneuploidies in the Danish iPSYCH2015 dataset: a case-cohort study.

BACKGROUND: Increased prevalence of mental illness has been reported in clinical studies of sex chromosome aneuploidies, but accurate population-based estimates of the prevalence and clinical detection rate of sex chromosome aneuploidies and the associated risks of psychiatric disorders are needed. In this study, we provide such estimates, valid for children and young adults of the contemporary Danish population. METHODS: We used the iPSYCH2015 case-cohort dataset, which is based on a source population of single-born individuals born in Denmark between May 1, 1981, and Dec 31, 2008. The case sample comprises all individuals from the source population with a diagnosis of any index psychiatric disorder (schizophrenia spectrum disorder, bipolar disorder, major depressive disorder, autism spectrum disorder, or ADHD) by the end of follow-up (Dec 31, 2015), registered in the hospital-based Danish Psychiatric Central Research Register. The cohort consists of individuals randomly selected from the source population, and overlaps with the case sample. Biobanked blood samples for individuals in the case and cohort samples underwent genotyping and quality-control filtering, after which we analysed microarray data to detect sex chromosome aneuploidy karyotypes (45,X, 47,XXX, 47,XXY, and 47,XYY). We estimated the population-valid prevalence of these karyotypes from the cohort sample. Weighted Cox proportional hazards models were used to estimate the risks of each index psychiatric disorder associated with each sex chromosome aneuploidy karyotype, by use of date of first hospitalisation with the index disorder in the respective case group and the cohort as outcome. The clinical detection rate was determined by comparing records of clinical diagnoses of genetic conditions from the Danish National Patient Register with sex chromosome aneuploidy karyotype determined by our study. FINDINGS: The assessed sample comprised 119 481 individuals (78 726 in the case sample and 43 326 in the cohort) who had genotyped and quality-control-filtered blood samples, including 64 533 (54%) people of gonadal male sex and 54 948 (46%) of gonadal female sex. Age during follow-up ranged from 0 to 34·7 years (mean 10·9 years [SD 3·5 years]). Information on ethnicity was not available. We identified 387 (0·3%) individuals as carriers of sex chromosome aneuploidies. The overall prevalence of sex chromosome aneuploidies was 1·5 per 1000 individuals. Each sex chromosome aneuploidy karyotype was associated with an increased risk of at least one index psychiatric disorder, with hazard ratios (HRs) of 2·20 (95% CI 1·42-3·39) for 47,XXY; 2·73 (1·25-6·00) for 47,XXX; 3·56 (1·01-12·53) for 45,X; and 4·30 (2·48-7·55) for 47,XYY. All karyotypes were associated with an increased risk of ADHD (HRs ranging from 1·99 [1·24-3·19] to 6·15 [1·63-23·19]), autism spectrum disorder (2·72 [1·72-4·32] to 8·45 [2·49-28·61]), and schizophrenia spectrum disorder (1·80 [1·15-2·80] to 4·60 [1·57-13·51]). Increased risk of major depressive disorder was found for individuals with 47,XXY (1·88 [1·07-3·33]) and 47,XYY (2·65 [1·12-5·90]), and of bipolar disorder for those with 47,XXX (4·32 [1·12-16·62]). The proportion of sex chromosome aneuploidy carriers who had been clinically diagnosed was 93% for 45,X, but lower for 47,XXY (22%), 47,XXX (15%), and 47,XYY (15%). Among carriers, the risk of diagnosis of at least one index psychiatric disorder did not significantly differ between those who had and had not been clinically diagnosed with sex chromosome aneuploidies (p=0·65). INTERPRETATION: Increased risks of psychiatric disorders associated with sex chromosome aneuploidies, combined with low rates of clinical diagnosis of sex chromosome aneuploidies, compromise the adequate provision of necessary health care and counselling to affected individuals and their families, which might be helped by increased application of genetic testing in clinical settings. FUNDING: Lundbeck Foundation and National Institutes of Health.

Accuracy of haplotype estimation and whole genome imputation affects complex trait analyses in complex biobanks.

Sample recruitment for research consortia, biobanks, and personal genomics companies span years, necessitating genotyping in batches, using different technologies. As marker content on genotyping arrays varies, integrating such datasets is non-trivial and its impact on haplotype estimation (phasing) and whole genome imputation, necessary steps for complex trait analysis, remains under-evaluated. Using the iPSYCH dataset, comprising 130,438 individuals, genotyped in two stages, on different arrays, we evaluated phasing and imputation performance across multiple phasing methods and data integration protocols. While phasing accuracy varied by choice of method and data integration protocol, imputation accuracy varied mostly between data integration protocols. We demonstrate an attenuation in imputation accuracy within samples of non-European origin, highlighting challenges to studying complex traits in diverse populations. Finally, imputation errors can bias association tests, reduce predictive utility of polygenic scores. Carefully optimized data integration strategies enhance accuracy and replicability of complex trait analyses in complex biobanks.

Interplay of ADHD Polygenic Liability With Birth-Related, Somatic, and Psychosocial Factors in ADHD: A Nationwide Study.

OBJECTIVE: Attention deficit hyperactivity disorder (ADHD) is a multifactorial neurodevelopmental disorder, yet the interplay between ADHD polygenic risk scores (PRSs) and other risk factors remains relatively unexplored. The authors investigated associations, confounding, and interactions of ADHD PRS with birth-related, somatic, and psychosocial factors previously associated with ADHD. METHODS: Participants included a random general population sample (N=21,578) and individuals diagnosed with ADHD (N=13,697) from the genotyped Danish iPSYCH2012 case cohort, born between 1981 and 2005. The authors derived ADHD PRSs and identified 24 factors previously associated with ADHD using national registers. Logistic regression was used to estimate associations of ADHD PRS with each risk factor in the general population. Cox models were used to evaluate confounding of risk factor associations with ADHD diagnosis by ADHD PRS and parental psychiatric history, and interactions between ADHD PRS and each risk factor. RESULTS: ADHD PRS was associated with 12 of 24 risk factors (odds ratio range, 1.03-1.30), namely, small gestational age, infections, traumatic brain injury, and most psychosocial risk factors. Nineteen risk factors were associated with ADHD diagnosis (odds ratio range, 1.20-3.68), and adjusting for ADHD PRS and parental psychiatric history led to only minor attenuations. Only the interaction between ADHD PRS and maternal autoimmune disease survived correction for multiple testing. CONCLUSIONS: Higher ADHD PRS in the general population is associated with small increases in risk for certain birth-related and somatic ADHD risk factors, and broadly to psychosocial adversity. Evidence of gene-environment interaction was limited, as was confounding by ADHD PRS and family psychiatric history on ADHD risk factor associations. This suggests that the majority of the investigated ADHD risk factors act largely independently of current ADHD PRS to increase risk of ADHD.

Cross-trait assortative mating is widespread and inflates genetic correlation estimates.

The observation of genetic correlations between disparate human traits has been interpreted as evidence of widespread pleiotropy. Here, we introduce cross-trait assortative mating (xAM) as an alternative explanation. We observe that xAM affects many phenotypes and that phenotypic cross-mate correlation estimates are strongly associated with genetic correlation estimates (R2=74%). We demonstrate that existing xAM plausibly accounts for substantial fractions of genetic correlation estimates and that previously reported genetic correlation estimates between some pairs of psychiatric disorders are congruent with xAM alone. Finally, we provide evidence for a history of xAM at the genetic level using cross-trait even/odd chromosome polygenic score correlations. Together, our results demonstrate that previous reports have likely overestimated the true genetic similarity between many phenotypes.

A comprehensive map of genetic relationships among diagnostic categories based on 48.6 million relative pairs from the Danish genealogy.

For more than half a century, Denmark has maintained population-wide demographic, health care, and socioeconomic registers that provide detailed information on the interaction between all residents and the extensive national social services system. We leverage this resource to reconstruct the genealogy of the entire nation based on all individuals legally residing in Denmark since 1968. We cross-reference 6,691,426 individuals with nationwide health care registers to estimate heritability and genetic correlations of 10 broad diagnostic categories involving all major organs and systems. Heritability estimates for mental disorders were consistently the highest across demographic cohorts (average h2 = 0.406, 95% CI = [0.403, 0.408]), whereas estimates for cancers were the lowest (average h2 = 0.130, 95% CI = [0.125, 0.134]). The average genetic correlation of each of the 10 diagnostic categories with the other nine was highest for gastrointestinal conditions (average rg = 0.567, 95% CI = [0.566, 0.567]) and lowest for urogenital conditions (average rg = 0.386, 95% CI = [0.385, 0.388]). Mental, pulmonary, gastrointestinal, and neurological conditions had similar genetic correlation profiles.

Associations between patterns in comorbid diagnostic trajectories of individuals with schizophrenia and etiological factors.

Schizophrenia is a heterogeneous disorder, exhibiting variability in presentation and outcomes that complicate treatment and recovery. To explore this heterogeneity, we leverage the comprehensive Danish health registries to conduct a prospective, longitudinal study from birth of 5432 individuals who would ultimately be diagnosed with schizophrenia, building individual trajectories that represent sequences of comorbid diagnoses, and describing patterns in the individual-level variability. We show that psychiatric comorbidity is prevalent among individuals with schizophrenia (82%) and multi-morbidity occur more frequently in specific, time-ordered pairs. Three latent factors capture 79% of variation in longitudinal comorbidity and broadly relate to the number of co-occurring diagnoses, the presence of child versus adult comorbidities and substance abuse. Clustering of the factor scores revealed five stable clusters of individuals, associated with specific risk factors and outcomes. The presentation and course of schizophrenia may be associated with heterogeneity in etiological factors including family history of mental disorders.

Genetic, Clinical, and Sociodemographic Factors Associated With Stimulant Treatment Outcomes in ADHD.

OBJECTIVE: Stimulant medications are effective for treating attention deficit hyperactivity disorder (ADHD), yet discontinuation and switch to nonstimulant ADHD medications are common. This study aimed to identify genetic, clinical, and sociodemographic factors influencing stimulant treatment initiation, discontinuation, and switch to nonstimulants in individuals with ADHD. METHODS: The authors obtained genetic and national register data for 9,133 individuals with ADHD from the Danish iPSYCH2012 sample and defined stimulant treatment initiation, discontinuation, and switch from prescriptions. For each stimulant treatment outcome, they examined associations with polygenic risk scores (PRSs) for psychiatric disorders and clinical and sociodemographic factors using survival analyses, and conducted genome-wide association studies (GWASs) and estimated single-nucleotide polymorphism heritability (h2SNP). RESULTS: Eighty-one percent of the sample initiated stimulant treatment. Within 2 years, 45% discontinued stimulants and 15% switched to nonstimulants. Bipolar disorder PRS (hazard ratio=1.05, 95% CI=1.02, 1.09) and schizophrenia PRS (hazard ratio=1.07, 95% CI=1.03, 1.11) were associated with discontinuation. Depression, bipolar disorder, and schizophrenia PRSs were marginally but not significantly associated with switch (hazard ratio range, 1.05-1.07). No associations were observed for ADHD and autism PRSs. Individuals diagnosed with ADHD at age 13 or older had higher rates of stimulant initiation, discontinuation, and switch (hazard ratio range, 1.27-2.01). Psychiatric comorbidities generally reduced rates of initiation (hazard ratio range, 0.84-0.88) and increased rates of discontinuation (hazard ratio range, 1.19-1.45) and switch (hazard ratio range, 1.40-2.08). h2SNP estimates were not significantly different from zero. No GWAS hits were identified for stimulant initiation or discontinuation. A locus on chromosome 16q23.3 reached genome-wide significance for switch. CONCLUSIONS: The study findings suggest that individuals with ADHD with higher polygenic liability for mood and/or psychotic disorders, delayed ADHD diagnosis, and psychiatric comorbidities have a higher risk for stimulant treatment discontinuation and switch to nonstimulants. Despite the study's limited sample size, one putative GWAS hit for switch was identified, illustrating the potential of utilizing genomics linked to prescription databases to advance ADHD pharmacogenomics.